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Abstract. Pseudorandom PSK pp enables parallel communication on 
the same carrier frequency and at the same time. We propose different 
signal processing methods to receive data modulated with pseudorandom 
PSK. This includes correlation with the carrier frequency which can be 
applied to signals in the kHz to MHz range and signal processing in the 
intermediate frequency where the correlation with the carrier frequency 
is performed analogous in the RF front end. We analyze the computation 
complexity for signal processing with the parameters of symbol length T 
and number of repetitions of each symbol K with pseudorandom PSK 
and show that the number of operations for each sampling point is 0 (K). 
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1 Introduction 


The proposed scheme for radio transmission of [T] repeats each symbol K times 
with pseudo random phase shifts at each repetition. If the receiver uses the same 
sequence of pseudo random phase shifts for correlation, it can demodulate the 
data from the signal. On the other hand, if a receiver uses another sequence of 
pseudo random phase shifts which is not applied by the sender, it receives noise 
with a smaller amplitude. This reduces the correlation of parallel and interfering 
transmissions on the same carrier frequency in expectation and makes possible 
to have parallel transmission at the same time, on the same carrier frequency, 
and possibly in the same spatial area. This is particularly interesting for ad hoc 
networks and distributed medium access control (MAC). 

According to [T], we can correlate the received signal y for known synchro¬ 
nization with a carrier signal with frequency / with 

T-l K-l 

p(t) = Y1 ■'^y{k-T + tk{t)+t + T) (1) 

t=0 k=0 

where 4>k are pseudorandom phase shifts with 0 < k < K, sample rate R, T 
samples per symbol, and 



if t < r - ^ 
otherwise 


tk (t) = 


( 2 ) 
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In a real channel, we also encounter effects like attenuation due to path path 
and multi-path propagation which is not subject of the analysis in this paper. 
Here, we show different signal processing methods how we can also efficiently 
synchronize to the transmitter signal besides demodulation of data after syn¬ 
chronization. This makes necessary correlate the received signal with the carrier 
(which is a sinusoidal carrier with pseudorandom phase shifts) at every sample 
position to determine the point of synchronization, i.e. start and end of a symbol. 

1.1 Contribution 

We propose correlation methods for the K-repetition scheme [T] necessary for 
synchronization and demodulating data. We improve the number computations 
of a correlation at a sample time from 0 {K ■ T) to 0 (T), where K is the num¬ 
ber of repetitions of each symbol and T is the number of samples per symbol. 
Table o shows the computational complexity of the proposed methods. We 
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phase error 

carrier 

P 
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0{T) 

0 

carrier 
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Table 1. Correlation methods with carrier frequency / with T samples respectively p 
periods per symbol and K repetitions of each symbol with pseudo-random phase shifts 


present methods for signal processing at the carrier frequency or at a smaller 
intermediate frequency, where correlation with the carrier frequency is done anal¬ 
ogously in the radio fronted. For sampling at the carrier frequency we assume 
oversampling with T samples per symbol and consider the correlation at each 
sample point. The phase error relates only to computational errors of the signal 
processing method and not e.g. errors due to sampling. 

2 Correlation with Carrier Signal 

We show algorithms for signal processing with correlation of the received signal 
y (r) at sampling point r and a reference signal containing a sinusoidal carrier 
with frequency /. We suppose oversampling with a rate R with T samples per 
symbol. 

In the first algorithm pmai {t) for correlation at sample point r, we set K = 1 
and get PSK without repetition of the symbols. We assume that at time t = 0, 
the phase angle of the reference signal for correlation is zero with = 0. 

We intend to compute the correlation of a window of T samples for each sample 
point, i.e. for the windows with sample ranges [0, T — 1], [1, T],..., [r, r -|- T — 1]. 
If we compare the correlation sum at two subsequent points r — 1 and r with win¬ 
dows [r — r, r — 1] and [r — T -|- 1, r], the calculation only differs by the sample 
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points t — T and r and the remaining T — 1 summands are the same. Comparable 
to a moving average (ma), we can compute the correlation for the received signal 
y (r) at sample point r from the correlation result Pmai (t — 1 ) at the preceding 
sample point with 

Pmai (r) = p„,ai {t - 1) + v (t) • - y {t - T) ■ e-^ 2 ./.(r-r)/fl _ ( 3 ) 

This equation is recursive and we can initialize Pmai (0) := 0 and y (r) — 0 for r < 
0. If we buffer the intermediate results of x (t) j-j^g buffer with T 

elements (which is initialized with zeroes), we can save the computational costs 
for the third term x (r) • q^ben, we only need one multiplication 

and two additions of complex values for each sample point. Since we correlate 
the received signal y (r) with a periodically carrier , the complex values 

can be pre-calculated and kept in a lookup table for given discrete sampling 
points. 

Lemma 1. After T sequential runs of pmai (t) at t = 0,... ,T—1, the operation 
Pmai outputs the Same result as p (r) of Equation 0 for K = 1. 

Proof: For K = 1, Equation 0 simplifies to 

T-l 

p (t) = ^ y (t + t — T + 1 ) • ^ (^ 4 ^ 

t=o 

The original equation Q returns at t = 0 which requires the reception of the 
samples at r = [0, T) for simplihcation. Without loss of generality we have added 
here the term — (T — 1) and get the first correlation output at r = T — 1. For 
the first T runs, the term y (r — T) of Equation ® is zero. With pmai (0) := 0, 
we get after T iterations a sum of y (t) • with r = [0, T) which is 

T-l 

Pmai (T - 1 ) = ^ y (t) . = p{T-l) . 

The proof for T > 0 follows by induction. 

Pmai T 1) = Pmai ('^) 

+y (r + 1 ) • _ y ((r + 1 ) - T) • 

T-l 

^ y (r + t - r + 1 ) • 

+y (r + 1) • e-J^^f<r+l)/R _ y + 4 ) _ . ^- 32 ^f-((r+l)-T)/R 

T 

= '^y{r + t-T + l)- 
— Pmai ('^ T 1 ) 
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For a repetition of symbols with iC > 1, we can use a moving average as well 
with 


K-l 


Pma2 


(t) = ^ Pmal {t -T ■ k) -e 


.j4>k 


( 5 ) 




For each summand in Equation ([^ we need two additions and two multipli¬ 
cations of complex values. The overall computational complexity of (/)am 2 {t) is 
0{K). 


Lemma 2. After T ■ K sequential runs of Pma 2 {T) t = 0,... ,T — 1, the 
operation Pma 2 outputs the same result as p ir) of Equation 0 for K > 1. 

Proof: We have already proven the statement for if = 1 in the proof of 
Lemma To proof the case K > 1 we insert Equation Q into ([^ and get 

K-l /T-l \ 

Pma2(r) = ^y(T-r.fc + t-T + l).e-^W/« • . 

k—0 \ t—0 / 


We state in [T] that the delay tk ft) replaces the phase shift with operation 
and hence Pma .2 (t) = p (r). Please note, that we assumed for the computation of 
p (t) to have the samples for [ t , t + K ■ T) in advance and in pnia 2 we compute 
the result from the preceding samples {t — K • T, r] which taken into account 
in the last transformation of the equation. Consequently, we need T ■ K itera¬ 
tions before Pma 2 (t) is initialized and outputs the correlation of a window of 
length T ■ K. □ 


In the next calculation procedure (fmaz (t), we reduce the number of multi¬ 
plications from 2 ■ K to one with the sacrifice of a precision loss. 


Pma3 (^) — Pma3 i) 


/K-l 


\fe =0 / 




( 6 ) 


with pre-calculated values 


j./ _ 4>k ' R 


p) 


We assume that we compute the correlation for each sampling point r = 0,1,... 
sequentially. The third summand of Equation § has already been computed 
for the correlation at (T — 1) sample points before and we can buffer these in¬ 
termediate results in a ring buffer with T elements to reduce the number of 
computational steps. Then we only have to compute the sum of the second 
summand at sample points r which are K additions, perform only one multipli¬ 
cation with for correlation (instead of 2 • iF multiplications before). 
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and compute two additional sums of complex values for the overall sum of the 
three summands of Equation ([^ . 

We include the phase shifts (j)k of the repeated symbols with the time shift 
in Equation Q . The difference to Equation ([^ is that we do not perform the 
modulo-operation and we will make a computational error at the end of each 
symbol where the the sample point r -|- t'f. is in the following sample (instead 
of at the beginning of the correct symbol with the modulo operation). Let us 
estimate this error and assume that each symbol is send for p periods of the 
carrier frequency, i.e. for time p/ f for carrier frequency /. In worst case, = 0 
and = 2tt for 2 < k' < K. Then we read almost all information of the last 
period from the subsequent next symbol. We can represent the correlation of 
p periods as superposition of the correlation of p individual correlations, where 
we know that the last is erroneous. Eigure[2 shows an example of a correlation 





error margin 


Fig. 1. Error margin of correlation method pma 3 in Equation § for P = 4 periods 


with p = 4 periods of the carrier where each vector illustrates the correlation of 
one period and the last correlation with the red-colored vector is erroneous. The 
error in the amplitude is at most l/(p — 2) if the erroneous correlation of the 

last period has a phase shift of tt. The phase error is at most sin ^ 

for p > 2. With the sacrifice of the computation error we can first compute the 
sum of K samples and then perform the multiplication with the complex value 
to compute the correlation. Compared to Equation ([^ , this reduces the number 
of multiplications of a correlation from 2 • fC to 1. 


3 Processing at an Intermediate Preqnency 

Eor carrier frequencies in the GHz range, oversampling and correlation with a 
rate larger than the carrier frequency is mostly too fast for current hardware 
and correlation is implemented in hardware. The output of the radio front end 
are complex values in an intermediate frequency which is slower than the carrier 
frequency. Let us assume that we have G complex-valued samples pc (r) in the 
intermediate frequency. The discrete sample points r are in this case in the 
intermediate frequency. 
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Just like in Equation ([^ for a correlation with the carrier, we can use once 
again a moving average at the intermediate frequency for K = 1 with 

Pma4 (r) = /9ma4 (t - 1) + J/c (f) - TJc {t - G) . (8) 

Since the correlation is performed in analog hardware in the radio front end, this 
only needs two (complex) additions. 

For K > 1 we can use K moving averages with 


K-l 


PmaS (f) = > ^ Pma4 ij - k ■ G) ■ & 


J<Pk 




( 9 ) 


A correlation with Pmas (t) needs 3 ■ K — 1 complex additions and K complex 
multiplications for the phase shifts with 


yc{T - K ■ G) yc(T - 2G) y^G - G) ycir) 



Fig. 2. Block diagram for signal processing of method pmas (r) at sample point r for 
K repetitions of each symbol with pseudorandom phase shifts (f>k with 0 < k < K 


The block diagram in Figure]^ shows the computation of PmaS {t). The com¬ 
putations in the grey boxes are independent from each other and can be com¬ 
puted in parallel. The total sum (at the bottom of the chart) depends on the 
intermediate results (in the grey boxes) and it is also possible to add these K 
values with a tree-like concatenation with only log (AT) instead of K computa¬ 
tional steps shown in the chart. Then the overall number of computational steps 
is only & (log AT) if we can perform K operations in parallel. 
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4 Conclusions 

We show that computational costs for synchronizing and demodulating data in 
the K-repetition scheme are only increased by factor 0 (K) compared to sig¬ 
nal processing of a modulation without repeating each symbol K times. This 
indicates the practical applicability of the K-repetition scheme. We show meth¬ 
ods for correlation with the carrier (a sinusoidal signal with carrier frequency 
/ and K pseudorandom phase shifts). These are mostly applicable for carrier 
frequencies in the kHz to MHz range where the clock speed of a processor is high 
enough. We also show methods for demodulating the K-repetition scheme at an 
intermediate frequency where correlation with the carrier is done at the radio 
front. This is the typical application for carrier frequencies in the GHz range 
and needs 0 (K) operations per sampling point. If we can run K operations in 
parallel, we also show that the number of computational steps is only 0 (logKT). 
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